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I. INTRODUCTION 



The exclusive two-body decays of heavy quarkonium into hght hadrons have been studied 
in the framework of perturbative QCD by many authors (for reviews, see These 
processes exhibit a large hierarchy between the heavy quark mass, which sets the scale for 
annihilation processes, and the scales that determine the dynamical structure of the particles 
in the initial and final states. The large energy released in the annihilation of the heavy 
quark- ant iquark pair and the kinematics of the decay — with the products flying away 
from the decay point in two back-to-back, almost light-like directions — allow for rigorously 
deriving a factorization formula for the decay rate at leading twist (for an up-to-date review 
of the theoretical and experimental status of the exclusive decays into light hadrons, see j^). 

For the bottomonium system, a particularly interesting class of two-body final states 
is the ones containing two charmed mesons. In these cases the picture is complicated by 
the appearance of an additional intermediate scale, the charm mass rric, which is much 
smaller than the bottom mass rrib but is large enough to be perturbative. These decays 
differ significantly from those involving only light quarks. The creation of mesons that are 
made up of purely light quarks involves creating two quark-antiquark pairs, with the energy 
shared between the quark and antiquark in each pair. In the production of two D mesons, 
however, almost all the energy of the bottomonium is carried away by the heavy c and c, 
while the light quark and antiquark, which bind to the c and c respectively, carry away 
(boosted) residual energies. 

The existence of well-separated scales in the system and the intuitive picture of the decay 
process suggest to tackle the problem using a sequence of effective field theories (EFTs) 
that are obtained by subsequently integrating out the dynamics relevant to the perturbative 
scales rrih and rric. 

In the first step, we integrate out the scale rrify by describing the b and b with Non- 
Relativistic QCD (NRQCD) jj], and the highly energetic c and c with two copies of Soft- 
Collinear Effective Theory (SCET) j?] in opposite light-cone directions. In 

the second step, we integrate out the dynamics manifested at scales of order rric by treat- 
ing the quarkonium with potential NRQCD (pNRQCD) jl^ 12 L and the D mesons 
with a boosted version of Heavy-Quark Effective Theory (HQET) 13| IJ] 15| 16| 17 1 
The detailed explanation of why the aforementioned EFTs are employed is of- 
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fered in Sec. [TTl We will prove that, at leading order in the EFT expansion, the decay 
rate factors into a convolution of two perturbative matching coefficients and three (one for 
each hadron) non-perturbative matrix elements. The non-perturbative matrix elements are 
process-independent and encode information on both the initial and final states. 

For simplicity, in this paper we focus on the decays of the C-even quarkonia Xbj and 
rjb that, at leading order in the strong coupling a^, proceed via the emission of two virtual 
gluons. The same method can be generalized to the decays of C-odd states T and hb, which 
require an additional virtual gluon. We also refrain from processes that have vanishing 
contributions at leading order in the EFT power counting. So the specific processes studied 
in this paper are Xbo,2 ~* DD, Xbo.2 D*D*, and r]f, DD* + c.c. However, the EFT 
approach developed in this paper enables one to systematically include power-suppressed 
effects, making it possible to go beyond the leading-twist approximation. 

The study of the inclusive and exclusive charm production in bottomonium decays and 
of the role pla yed by the charm mass ttlc in such processes have recently drawn renewed 
attention [20| 2l| 22| 23|, in connection with the experimental advances spurred in the 
past few years by the abundance of bottomonium data produced at facilities like BABAR, 
BELLE, and CLEO. The most notable result was the observation of the bottomonium 

Furthermore, the 



ground state rjb-, recently reported by the BABAR collaboration 
CLEO collaboration published the first results for several exclusive decays of Xb into light 
hadrons and for the inclusive decay of Xb into open charm 26 1. In particular, they 
measured the branching ratio B{xbj where J is the total angular momentum 

of the Xb state, and conclusively showed that for J = 1 the production of open charm is 
substantial: B{xbii^P) D^X) = 12.59 ± 1.94%. For the J = 0,2 states the data are 
weaker, but the production of open charm still appears to be relevant. The measurements 
of the CLEO collaboration are in good agreement with the prediction of Bodwin et al. j2o| . 
where EFT techniques (in particular NRQCD) were for the first time applied to study the 
production of charm in bottomonium decays. 

The double-charm decay channels analyzed here have not yet been observed, so one of 
our aims is to see if they may be observable given the current data. Unfortunately, the poor 
knowledge of the D-meson matrix elements prevents us from providing definitive predictions 
for the decay rates T{xbj DD), T{xbj D*D*), and T{r]b DD* + c.c). As we will 
show, these rates are indeed strongly dependent on the parameters of the D- and D*-meson 
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distribution amplitudes, in particular on their first inverse moments Ad and Ad*: the rates 
vary by an order of magnitude in the accepted ranges for Ad and Ad*- On the other hand, 
the factorization formula implies that these channels, if measured with sufficient accuracy, 
could constrain the form of the D-meson distribution amplitude and the value of its first 
inverse moment. In turn, the details of the D-meson structure are relevant to other D-meson 
observables, which are crucial for a model-independent determination of the CKM matrix 



elements \Vcd\ and \ Vc 



27|. 



This paper is organized as follows. In Sec. [TTl we discuss the degrees of freedom and the 
EFTs we use. In Sec. IIII Al we match QCD onto NRQCD and SCET at the scale 2m;,. The 
renormalization-group equation (RGE) for the matching coefficient is derived and solved in 
Sec. IIIIBI In Sec. IIV Al the scale rric is integrated out by matching NRQCD and SCET onto 
pNRQCD and bHQET. The renormalization of the low-energy EFT operators is performed 
in Sec. IIVB| with some technical details left to App. |Al The decay rates are calculated in 
Sec. |V] using two model distribution amplitudes. In Sec. |VT]we draw our conclusions. 



II. DEGREES OF FREEDOM AND THE EFFECTIVE FIELD THEORIES 

Several well-separated scales are involved in the decays of the C-even bottomonia 77^ and 
Xbj into two D mesons, making them ideal processes for the application of EFT techniques. 
The distinctive structures of the bottomonium (a heavy quark-antiquark pair) and the D 
meson (a bound state of a heavy quark and a light quark) suggest that one needs different 
EFTs to describe the initial and final states. 

We first look at the initial state. The rjb is the ground state of the bottomonium system. 
It is a pseudoscalar particle, with spin S = 0, orbital angular momentum L = 0, and total 
angular momentum J = 0. In what follows we will often use the spectroscopic notation 
^^^^Lj, in which the rjb is denoted by ^Sq. The XbJ is a triplet of states with quantum 
numbers ^Pj. The rji, and Xbj are non-relativistic bound states of a 6 quark and a b antiquark. 
The scales in the system are the b quark mass m^, the relative momentum of the bb pair rrihW, 
the binding energy nibw'^, and Aqcd, the scale where QCD becomes strongly coupled, w is 
the relative velocity of the quark-antiquark pair in the meson, and from the bottomonium 
spectrum it can be inferred that w"^ ~ 0.1. Since rrib ^ Aqcd, can be integrated out 
in perturbation theory and the bottomonium can be described in NRQCD. The degrees 
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of freedom of NRQCD are non-relativistic heavy quarks and antiquarks, with energy and 
momentum [E, \p\) of order {rribw'^ , rribw) , hght quarks and gluons. In NRQCD, the gluons 
can be soft {rribW , nibw) , potential {rrihw'^ , mj,w) , and ultrasoft (usoft) {mbw"^ , mt,w'^) . The 
NRQCD Lagrangian is constructed as a systematic expansion in l/m^, whose first few terms 
are 

/:nrqcd = ^M^i^o + ^ + ^ + ...U + xM^/^o-^-^ + ... X, 

\ Zj / / fc^ 1 1 Off I \ 1 1 Off ^11 Off I 

where ip and annihilate a b quark and a b antiquark respectively, and ■ ■ ■ denotes higher- 
order contributions in l/m?,. In NRQCD several mass scales are still dynamical and different 
assumptions on the hierarchy of these scales may lead to different power countings for 
operators of higher dimensionality. However, as long as w <^ 1, higher-dimension operators 
are suppressed by powers of w (for a critical discussion on the different power countings we 
refer to Q). 

NRQCD still contains interactions that can excite the heavy quarkonium far from its 
mass shell, for example, through the interaction of a non-relativistic quark with a soft gluon. 
In the case m^w ^ Aqcd, we can integrate out these fluctuations, matching perturbatively 
NRQCD onto a low-energy effective theory, pNRQCD. We are then left with a theory of non- 
relativistic quarks and ultrasoft gluons, with non-local potentials induced by the integration 
over soft- and potential-gluon modes. The interactions of the heavy quark with ultrasoft 
gluons are still described by the NRQCD Lagrangian, with the constraint that all the gluons 
are ultrasoft. In the weak coupling regime rrihW ^ Aqcd, the potentials are organized by an 
expansion in as{miyw), l/m^, and r, where r is the distance between the quark and antiquark 
in the quarkonium, r ~ l/rrihW. If we assume mhw"^ ~ Aqcd, each term in the expansion 
has a definite power counting in w and the leading potential is Coulombic V ~ as{mbw)/r. 

An alternative approach, which does not require a two-step matching, has been developed 



in the effective theory vNRQCD |28| [29 1 



301 ] 31| . In the vNRQCD approach there is only 



one EFT below rrib, which is obtained by integrating out all the off-shell fluctuations at the 
hard scale rrih and introducing different fields for various propagating degrees of freedom 
(non-relativistic quarks and soft and ultrasoft gluons). In spite of the differences between 
the two formalisms, pNRQCD and vNRQCD give equivalent final answers in all the known 
examples in which both theories can be applied. 

We now turn to the structure of the D meson. The most relevant features of the D meson 
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are captured by a description in HQET. In HQET, in order to inte 



rric, the momentum of the heavy quark is generically written as 



grate out the inert scale 



15| 



p = rricV + , (1) 

where v is the four-velocity label, satisfying v'^ = 1, and k is the residual momentum. If 
one chooses v to be the center-of-mass velocity of the D meson, k scales as ~ f Aqcd- 
Introducing the light-cone vectors n^^ = (1,0,0,1) and = (1,0,0,-1), one can express 
the residual momentum in light-cone coordinates, k^ = n - kn^/2 + n-kn^ /2 + k^^ or simply 
k = {n-k,n-k, k^). There are two relevant frames. One is the D-meson rest frame, in which 

V is conveniently chosen as vq = (1,0,0,0), and the other is the bottomonium rest frame, 
in which the D mesons are highly boosted in opposite directions, with v chosen as v = ve>, 
the four- velocity of one of the D mesons. By a simple consideration of kinematics and the 
scaling k ~ f Aqcd, one can work out the scalings for k in the two frames. In the D-meson 
rest frame, k ~ Aqcd(15 1, 1), and in the bottomonium rest frame (supposing the D meson 
moving in the positive 2;-direction) , 

k ~ Aqcd {n ■VD.fi- vd, 1) ~ Aqcd^ ■ vd (A^, 1, A) , (2) 

where n ■ vd ^ 2mb/mc and A = mc/2mh <^ 1. It is convenient for the calculation in this 
paper to use the bottomonium rest frame, so we drop the subscript in vd and we assume 

V = Vd in the rest of this paper. The momentum scaling in Eq. ([2]) is called ultracollinear 
(ucollinear), and boosted HQET (bHQET) is the theory that describes heavy quarks with 
ultracollinear residual momenta and light degrees of freedom (including gluons and light 
quarks) with the same momentum scaling. 

The bHQET Lagrangian is organized as a series in powers of Aqcd/'^c and, for residual 



momentum ultracollinear in the n-direction, the leading term is 
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-CbHQET = hniv ■ Dhn , (3) 

where the field /i„ annihilates a heavy quark and the covariant derivative D contains ultra- 
collinear and ultrasoft gluons, 

iD^ = — {in ■ d + gn ■ An) + — {in ■ d + gn ■ + gn ■ A^s) + [id^^ + . (4) 

The ultrasoft gluons only enter in the small component of the covariant derivative. This 
fact can be exploited to decouple ultrasoft and ultracollinear modes in the leading-order 
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TABLE I: Degrees of freedom in EFTi(NRQCD + SCET). w is the bb relative velocity in the 
bottomonium rest frame, while A ~ mc/2mb is the SCET expansion parameter. We assume 
rubW ~ rric (or, equivalently, ~ A) and rribVLp' ~ nibX^ ~ Aqcd- 



Lagrangian through a field redefinition reminiscent of the colhnear-ultrasoft decouphng in 
SCET [Tf ISj. The ultracoUinear-ultrasoft decouphng is an essential ingredient for the 



factorization of the decay rate. 

Therefore, the appropriate EFT to calculate the decay rate is a combination of pNRQCD, 
for the bottomonium, and two copies of bHQET, with fields collinear to the n and n direc- 
tions, for the D and D mesons, symbolically written as EFTn = pNRQCD + bHQET. 

As we mentioned earlier, we plan to describe the bottomonium structure with a two-step 
scheme QCD — NRQCD pNRQCD. However, at the intermediate stage, where we first 
integrate out the hard scale 2mb and arrive at the scale m^w, the D meson cannot yet be 
described in bHQET. This is because the interactions relevant at the intermediate scale mbW 
can change the c-quark velocity and leave the D meson off-shell of order ~ {mbwY ~ ^1 S> 
Aqqj-,. Highly energetic c and c travelling in opposite directions can be described properly 
by SCET with mass. Thus, at the scale = 2mfe, we match QCD onto an intermediate 
EFT, EFTi = NRQCD + SCET, in which the EFT expansion is organized by A and w. The 
degrees of freedom of EFTi are tabulated in Tab. [B 

Then, we integrate out irtc and nibW at the same time, matching EFTj onto EFTn at the 
scale /i' = nic. In EFTn, the low-energy approximation is organized by Aqcd/'^c and w. 
The degrees of freedom of EFTn are summarized in Tab. [Tll When no subscript is specified 
in the rest of this paper, any reference to EFT applies to both EFTi and EFTn. To facilitate 
the power counting, we adopt w ~ A ~ Aqcd/''^c- As a first study, we will perform in this 
paper the leading-order calculation of the bottomonium decay rates. 
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TABLE II: Degrees of freedom in EFTii(pNRQCD + bHQET). The scale Q in bHQET is Q = 
n ■ w'Aqcd for the n-collinear sector and Q = n ■ vAqqd for the n-collinear sector, n ■ v' and 
n ■ V are the large light-cone components of the D-meson velocities in the bottomonium rest frame, 
n ■ v' ^ n ■ V ^ 2nii,/mc- A and w are defined as in Tab. HI 




FIG. 1: Matching QCD onto EFTj. On the r.h.s., the double lines represent the non-relativistic b 
(6) (anti)quark, while the dashed lines represent the collinear c (c) (anti)quark. 

III. NRQCD + SCET 
A. Matching 

In the first step, we integrate out the dynamics related to the hard scale 2mb by matching 
the QCD diagrams for the production of a c c pair in the annihilation of a 6 6 pair onto their 
EFTi counterparts. The tree-level diagrams for the process are shown in Fig. [H The gluon 
propagator in the QCD diagram has off-shellness of order = (2mfe)^ and it is not resolved 
in EFTi, giving rise to a point-like interaction. 

We calculate the diagrams on shell, finding 

iJqcD = «C(/u) Jefti(/u) , (5) 

with, at tree level, 

JEFT, = xyjybrnSll,±t^S^X^ and = 2m,) = ^^^^^ , (6) 
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where t"' are color matrices and the symbol o"'^ denotes the four matrices cr^ = (!,<?), with 
a the Pauli matrices. The subscript ± refers to the components orthogonal to the light- 
cone vectors and n^. The fields ipb and xl ^i^re two-component spinors that annihilate 
respectively a b quark and a b antiquark. Xn,n-p Xn n-p collinear gauge-invariant 
fermion fields: 

Xn,n-p = ^yWn^n)n-p i Xn,n-p = {^n^n)n-p ; (7) 

where Wn is defined as 

perms 

ly^^ has an analogous definition with n —>■ n. Collinear fields are labelled by the large 
component of their momentum. Note, however, we omit in Eq. ([6]) the subscripts n ■ p and 
n ■ p of the collinear fermion fields, in order to simplify the notation. The operator n ■ V in 
the definition ([8]) is a label operator that extracts the large component of the momentum of 
a collinear field, n ■ V (pn,n-p = n ■ p(pn,n-p, where 4>n,n-p is a generic collinear field. Sn(n) is a 
soft Wilson line, 

perms 

where the operator n ■ V acts on soft fields, n ■ V (ps = n ■ k (ps- 

Since in SCET different gluon modes are represented by different fields, we have to 
guarantee the gauge invariance of the operator Jefti under separate soft and collinear gauge 
transformations. A soft transformation is defined by Vs{x) = exp {ijSgt"-), with d^V ~ 
2mf,(A, A, A), while a gauge transformation U{x) is n-collinear if U{x) = exp {ia"'{x)t°') and 
d^U{x) ~ 2m{,(A^, 1, A). It has been shown in Ref. |7] that collinear fields do not transform 
under a soft transformation and that the combination PV^^n is gauge invariant under a 
collinear transformation. Soft fields do not transform under collinear transformations but 
they do under soft transformations. For example, the NRQCD quark and antiquark fields 
transform as ipb Vs{x)ip}j. The soft Wilson line has the same transformation, Sn — >■ 
Vs{x)Sn- Therefore, xl(^±'t°"^b transforms as an octet under soft gauge transformations. 
Since Xn ^nlfj.i-^"' ^nXn behaves like an octet as well, JefTi is invariant. It is worth noting 
that the soft Wilson lines are necessary to guarantee the gauge invariance of Jefti- We 
have explicitly checked their appearance at one gluon by matching QCD diagrams like the 
one in Fig. [H with all the possible attachments of an extra soft or collinear gluon, onto 
four-fermion operators in EFTi. 



B. Running 



The matching coefficient C and the effective operator JefTi depend on the renormahzation 
scale n- Since the effective operator is sensitive to the low-energy scales in EFTj, logarithms 
that would appear in the evaluation of Jefti are minimized by the choice /i ~ rric. On 
the other hand, since the coefficient encodes the high-energy dynamics of the scale 2mb, 
such a choice would induce large logarithms of mc/2m{, in the matching coefficient. These 
logarithms can be resummed using RGEs in NRQCD + SCET. 



The fi dependence of JefTj is governed by an equation of the following form [32 1 

d 



dlnfi 



-^EFTi(/i) = -7EFTi(/^)</eFTi(/^) 



(10) 



where the anomalous dimension 7efti is given by 



7efti 



d 



EFT 



^ dlnfi 



^EFTt 



(11) 



and 2'efti is the counterterm that relates the bare operator J^pT^ to the renormalized one, 
■^EFTi — ^EFTiifJ') JEFTiifJ') ■ Siucc the l.h.s. of Eq. (I5l) is independent of the scale /i, the 
RGE ffTOl) can be recast as an equation for the matching coefficient C(/i), 

d 



dlnfi 



C{fi) = 7EFTi(/i)<^(^) • 



(12) 



The counterterm Zeftj cancels the divergences that appear in Green functions with the 
insertion of the operator Jefti- We calculate 2'efti in the MS scheme by evaluating the 
divergent part of the four-point Green function at one loop, given by the diagrams in Figs. 

El-H 







FIG. 2: Soft diagrams at one loop. 



Since in NRQCD we do not introduce different gluon fields for different momentum modes, 
"soft" and "ultrasoft" in Fig. [2] and Fig. [3] refer to the convention that we impose soft or 
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FIG. 3: Ultrasoft diagrams at one loop. 





FIG. 4: Collinear diagrams at one loop. 

ultrasoft scaling to the corresponding loop momentum. The potential region, which should 
be considered in the diagrams of Fig. [2], does not give any divergent contribution. 

The integrals are evaluated in dimensional regularization, with d = 4 — 2e. We regulate 
the infrared divergences by keeping the non-relativistic b and b and the collinear c and c 



off-shell: Ef,^i - p^-j2mb = Ab, p\ 



and p\ 



m:. 



A^. We power count the 



c-quark off-shellness as A^ ~ A^ ~ m^A^ and the 6-quark off-shellness as A 



nibw'^. We 



also assume A^, A^ > 0. To avoid double counting, we define the one- loop integrals with 



the 0-bin subtraction 



33|. 



Even with an off-shellness, the soft diagrams in Fig. [2] do not contain any scale and they 
are completely cancelled by their 0-bin. 

The divergent part of the ultrasoft diagrams in Fig. [3] is 



iM 



usoft 



— - - In 

e-^ e 



A^A^ 



+ ^lln(-l-iO)- — i 



Jl 



EFT, 



(13) 



where Cp = (N^ — l)/2Ai'c and /i is the MS unit mass, /i^ = 47r/i^gexp {—'Je)- The first 
term in the curly brackets of Eq. (fT3|) corresponds to the sum of the divergences in the 
second diagram in Fig. [3l where an ultrasoft gluon is exchanged between the c and c quarks 
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collinear in back-to-back directions, and those in the last four diagrams of the same figure, 

which contain ultrasoft interactions between the initial and final states. The second term is 

an extra imaginary piece generated by the second diagram in Fig. [31 The —iO prescription 

in the argument of the logarithm, where is a positive infinitesimal quantity, follows from 

the prescriptions in the quark propagators and from the choice A^, > 0. The divergences 

arising from the ultrasoft exchanges between the b b pair in the first diagram in Fig. [3] are 

encoded in the last term in Eq. f|T3l) . 

The initial and final states cannot interact by exchanging collinear gluons because the 

emission or absorption of a collinear gluon would give the b quark an off-shellness of order 

ml, which cannot appear in the effective theory. For the same reason, the c and c cannot 

exchange n or n-coUinear gluons. The only collinear loop diagrams consist of the emission 

of a r;,(r;,)-collinear gluon from the Wilson line Wn{n) in </efTi and its absorption by the c(c) 

quark, as shown in Fig. HI The divergent part of the sum of the two collinear diagrams is 

'A^A^M 
In ( 1 1 Jefti . (14) 



iM coll = i—2CF 



4 + i|2 



The collinear diagrams are calculated with a 0-bin subtraction 33j, that is, we subtract from 
the naive collinear integrals the same integrals in the limit in which the loop momentum is 
ultrasoft. In this way we avoid double counting between the diagrams in Figs. [3] and HI 
Summing Eqs. ( |T3l) and ( |T4l) and adding factors of for each field. 



1 a. 

£271-' 



the divergent piece becomes 

. Oi 



iM 



div 



2 2 
+ - 



In 



e An 



1 ,r "^TT 1 

+ -N, + - — 

e e Nr. 



EFT, 



The counterterm Z-^fti is chosen so as to cancel the divergence in Eq. f|T5|) . 



^EFTt 



as 
Att 



2 2 /3 , 
— + - I - - In 



n ■ pcTi ■ pc 



1 ^ in 1 
-Nr + - — 

e € Nr 



(15) 



(16) 



From the definition ffTTl) . Eq. f|T6l) . and recalling that das/ din fi = —2eas + 0{al), the 
anomalous dimension at one loop is 



7efTi 



.2^iM<!3C^ + iV, + 4C^ln 



+^4 



^ , ... . (17) 

An y ^ ^ ^ \ y/n ■ pcU ■ Pc,' '^^ ' 

An important feature of the anomalous dimension (ITTI) is the presence of a term propor- 
tional to In/i. Because of this term, the RGE (fT2|) can be used to resum Sudakov double 
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logarithms. As we will show shortly, the general solution of Eq. (fT2|) can be written in the 
following form: 

= C{^l,) { ^ expf/(/io,/x) , (18) 

where g and U depend on the initial scale and the final scale /i that we run down to. For 
an anomalous dimension of the form (fT71) . U can be expanded as a series, 

oo n+1 



f/(/xo, /x) = V a:(/io) V Un,L ln"-^+' - • (19) 



ji=l L=0 

If /i//xo ^ 1; the most relevant terms in the expansion (fT9|) are those with L = 0, which we 
call "leading logs" (LL). Terms with higher L are subleading; we call the terms with L = 1 
"next-to- leading logs" (NLL), those with L = 2 "next-to-next-leading logs" (NNLL), and, 
if L = m, we denote them with N'^LL. The RGE (1121) determines the coefficients in the 
expansion f|T9|) . With the anomalous dimensions written as 

= -2 {■'(-■^ + lu {^^^^^^ } . (20) 

where 7(as) and T{as) are series in powers of as, 

it can be proved that the coefficients of the LL, Uno, are determined by the knowledge of F*-"^ 
and of the QCD /3 function at one loop. The NLL coefficients Uni are instead completely 
determined if F and (3 are known at two loops and 7(as) at one loop. 

In the case we are studying, the ratio of the scales fi/fiQ ^ mc/2mf, is not extremely 
small. Indeed, as to be seen shortly, the numerical contributions of the LL and NLL terms 
in the series f ll9l) are of the same size. It is therefore important to work at NLL accuracy, 
which requires the calculation of the coefficient of In /i to two loops. The factors of In /x are 



induced by cusp ang. 



es involving light-like Wilson lines and their coefficients are universal 



34| . The cusp anomalous dimension Fcusp(«s) is known at two loops j34| . 



rcusp(as) - . r^^^sp + ( . ) r^ulp > (21) 



47r ^"^P V47r7 

with 



p(0) _ Ary p(l) _ 

cusp ^'~'J^ , J- cusp ^'-V^ 



67 7r^\ 10 

] Uf 

9 3 / 9 ^ 



(22) 
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while the constant of proportionahty between T{as) and rcusp(«s) is fixed by the one-loop 
calculation. Since we have determined 7^°-*, 



7 



7r 



(23) 



and the (3 function is known, we have all the ingredients to provide the NLL approximation 
for U{fio, fi) and g{fio, /^)- Taking into account the tree-level initial condition in Eq. ([6]), Eq. 
(fT8|) determines the leading-order matching coefficient, with NLL resummation. 
The solution f|T8l) can be derived by writing Eq. f|T2l) as 



d\nC 



da 



7(a) + r 



cusp 



a 



In 



+ 



da' 



(24) 



where we have used the definition of the (3 function, l3{a) = da/dlnfi, to write In/i and 
dlnfi in terms of a. Integrating both sides from fiQ to fi and exponentiating the result we 
find the form given in Eq. (fTSl) . with 

""^^^^ da 



o r^^^ da ^ , ^ , , r da' ] 
U{fio,fi) = -2 -— <^7(a) + rcusp(a) / -^r^ \ 



^(/io,/i) = -2 / -^^rcusp(a) 



(25) 



At NLL, we find 



2ttT 



(0) 
cusp 



1-rlnr /?o7r, 



(0) 



+ 



27rr 



(0) 



■ Inr -|- 



In r 



+ 



cusp 
(0) 



.(1) o 
cusp yl 

Jo 



.(0) 
cusp 



1 — r -|- In r 
Air 



(26) 



/3o 



and 



-(0) 
- cusp 



(0) 
cusp 



Air 



(27) 



where r = as{fi) / asifJ^b) and we have renamed the initial scale fib, to denote its connection 
to the scale 2m;,. In Eqs. and (B7|l we have used the two-loop beta function. 



with 



p{as 



Po = ll - -rif , 



(28) 



(29) 



In Eq. (126|) we have kept the contributions of the real and imaginary part of 7'''''* separated. 
The imaginary part of 7'^'^^ changes the phase of the matching coefficient C{fi), but this phase 
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is irrelevant for the calculation of physical observables like the decay rate, which depend on 
the square modulus of C{fi). In Sec. |V]the factor U{jj,h^lj) will be evaluated between the 
scales /ift = 2mf, and /i = m^, with = 4 active quark flavors. The numerical evaluation 
shows that the LL term, represented by the first term in the brackets in Eq. fl2Bl) . is slightly 
smaller than and have the opposite sign of the term proportional to ^^^^ which dominates 
the NLL contribution. This observation confirms, a posteriori, the necessity to work at NLL 
accuracy in the resummation of logarithms of mc/2mh. 

The RGE ( IT2|) and its solution (fTS!) thus allow us to rewrite Eq. ([5]) as 

-^QCD = C'(/i)JEFTi(/i) = C{fib = 2mfe) expf/(2m6,mc)JEFTi(/i = rUc) , 

which avoids the occurrence of any large logarithm in the matching coefficient or in the 
matrix element of the effective operator. 

IV. pNRQCD + bHQET 
A. Matching 

In the second step, we integrate out the soft modes by matching EFTj onto EFTn. In 
NRQCD + SCET, contributions to the exclusive decay processes are obtained by considering 
time-ordered products of Jefti and the terms in the EFTi Lagrangian that contain soft- 
gluon emissions. The soft gluons have enough virtuality to produce a pair of light quarks 
travelling in opposite directions with ultracoUinear momentum scaling. These light quarks 
bind to the charm quarks to form back-to-back D mesons. The total momentum of two back- 
to-back ultracoUinear quarks is 2mfeAQCD/'"^c (1, 1, A) and the invariant mass of the pair is 

~ {2mbAQCD / fncY ~ ^l- in NRQCD + SCET, only soft gluons have enough energy to 
produce them. The time-ordered products in NRQCD + SCET are matched onto six-fermion 
operators in pNRQCD + bHQET, where fluctuations of order cannot be resolved. 

We consider the scale /i' = nic to be much bigger than Aqcd, so the matching can be done 
in perturbation theory. The Feynman diagrams contributing to the matching are shown in 
Fig. [51 The gluon and the 6-quark propagators have off-shellness of order ml, so the two 
diagrams on the l.h.s. match onto six-fermion operators on the r.h.s. 

The amplitude for the decay of a bottomonium with quantum numbers ^^^^Lj into two 
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FIG. 5: Matching NRQCD + SCET onto pNRQCD + bHQET. On the r.h.s. the double sohd hnes 
represent heavy b (b) (anti)quarks, the double dashed lines bHQET c (c) (anti)quarks, and the 
single dashed lines collinear light quarks. 

D mesons has the following form: 

J LO U 

(30) 

A and B, which label the final states and the EFTn operators denote the possible 

parity, spin, and polarization of the D mesons, A,B — {P, Vl, Vt}, indicating respectively 
a pseudoscalar D meson, a longitudinally-polarized vector meson D*, and a transversely- 
polarized vector meson D*. Unhke JefTi, we have dropped the subscript EFTn in O^^^^^-' 
in order to simplify the notation. 

The EFTn operators that contribute to the decay of the P-wave states are 



Lf,r-u ~ -i--r u ■ -n 2^ I ~ \ — ■- ■ I /v-n ^v-n- v" " ' ^ 2 



(31) 



where Ph]_cr'2 is a symmetric, traceless tensor. 



At leading order in the EFTn expansion, the rjb can only decay into a pseudoscalar and a 

vector meson, with an operator given by 



(32) 
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For later convenience, in the definition of tlie effective operators (13T!) and (132|) we have 
factored out the term which is related to the D-meson decay constant. The definition 

of will become clear when we introduce the D-meson distribution amplitudes. The 

fields Xn ^"^^ Xn ultracollinear gauge-invariant light-quark fields, while = W^h'^ and 
~ ^n^n bHQET hcavy-quark fields, which are invariant under an ultracollinear 
gauge transformation. The Wilson lines Wn and Wn have the same definition as in Eq. ([8]), 
with the restriction to ultracollinear gluons. Eqs. fl3Tl) and fl32|) allow us to interpret u as the 
component of the light-quark momentum along the direction n. Similarly, u represents the 
component of the light-antiquark momentum along n. The minus sign in the delta function 
6{—Lj — n -V) is chosen so that u is positive. 
The tree-level matching coefficients are 

T{u,u,iJ,iJ = nic; Pj) - 



iV2 mb u + uj 

T = m^; So) = — — — . 

mh 2uj + oj 

Note that, at leading order in the EFTn expansion, the matching coefficient 
T{u,uj, fi, fi'] •^Pj) is independent of the spin and polarization of the final states, or of the 
total angular momentum J of the Xb- 

An important feature of bHQET is that the ultracollinear and ultrasoft sectors can be 
decoupled at leading order in the power counting by a field redefinition reminiscent of the 
coUinear-usoft decoupling in SCET 0] [ist. For bHQET in the n direction, the decoupling 
is achieved by defining h'^ ^nh'^ and ilXh "where y„ is an ultrasoft Wilson line, 

y.= 5^ [^^v{-^n.A^s)] . (34) 

perms 

An analogous redefinition with n —>■ n decouples ultrasoft from fi-ultracollinear quarks and 
gluons. These redefinitions do not affect the operators in Eqs. (l3T!l and (132|) because all 
the induced Wilson lines cancel out. As a consequence, at leading order in the EFTn power 
counting, there is no interaction between the initial and the final states, since the former 
can only emit and absorb ultrasoft gluons that do not couple to ultracollinear degrees of 
freedom. Furthermore, fields in the two copies of bHQET, boosted in opposite directions, 
cannot interact with each other because the interaction with a n-ultracoUinear gluon would 
give a n-ultracoUinear quark or gluon a virtuality of order m^, which, however, cannot appear 
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in EFTn- The matrix elements of the operators Oj^^ {uj,uj, fx), therefore, factorize as 

^, f^'m ={0\xl TlT'^'Albb) {A\ni 5{-u-n.V) |0) 
{B\tj{uJ-n.V^) |rBK|0), 

where Va = {75, and = {1, pf, - a±, Pi±cr±}- The charge-conjugated contribu- 

tion is understood in the rjb case. 

The quarkonium state and the D mesons in Eq. fl35l) have respectively non-relativistic 
and HQET normalization: 

{Xu{E',p')\xtAE,p)) = {2iTf5^'\p-p') , {D{v',k')\D{v,k)) = 2vXy{27if6^'\k -k') , 

where f° is the 0th component of the 4- velocity v^. 

The D-meson matrix elements can be expressed in terms of the D-meson light-cone 
distribution amplitudes: 

{P\x^i7'6{uj-n.V^)nm =^Fp(/i')^0p(^,/x'), (36) 

J _ ft ' V 

{VL\xLl5{^-n-V^)nm =Fy,{^^')^<PyAu;,^i'), (37) 
{VT\xLillH^-^-^^)K\0) =FyM^e^^<l>y,iuj,fx') , (38) 

where is the transverse polarization of the vector meson. The constants FAifJ^'), with 
A = {P,Vl,Vt}, are related to the matrix elements of the local heavy-light currents in 
coordinate space. In the heavy-quark limit, where D and D* are degenerate. Fa is the same 
for all the three states: F = Fp = Fy^ = Fy^. In this limit, 

ml^l'KmP) = -^Fi^^')^ . (39) 

At tree level, the matrix element is proportional to the D-meson decay constant fp, = 
205.8 ± 8.5 ± 2.5 MeV |35|. More precisely, F{^') = foy/mD, where the factor y/mo is due 
to HQET normalization. The scale dependence of F is determined by the renormalization 
of heavy-light HQET currents. At one loop, Ref. |32] showed that 

^FM = -7.fM=3C,^FM. (40) 

The pNRQCD matrix elements can be expressed in terms of the heavy quarkonium 
wavefunctions. The operator xlpb ■ (^±'ipb contains a component with J = and a component 
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with J = 2 and Jz = 0, so its matrix element has non- vanishing overlap with both xw 
and Xb2- The operator xlPb^(^±'^b instead has only contributions with J = 2 and Jz = ±2 
and therefore it only overlaps with Xb2- In terms of the bottomonium wavefunctions, the 
pNRQCD matrix elements are expressed as 



{o\xlpb ■ ^±i'b\xbo) 
{^IxiPb ■ ff±A\xb2) 



v/3 V 27r ^-'"^ ' ' 



(£(2)+^(-2)^ 



2-K 



(41) 
(42) 
(43) 



where i?^^^(0) is the derivative of the radial wavefunction of the XbJ evaluated at the origin. 
At leading order, the pNRQCD Hamiltonian does not depend on J, so, up to corrections of 
order w^, R'^^^^i^) = R'xto^^^' '^^^ numerical pre-factors in Eqs. fj4Tl) and fj42|) follow from 
decomposing pb ■ a± into components with definite Jz- s^^l is the polarization tensor of the 
Xb2 state, and Eq. fH5]) states that, at leading order in the w'^ expansion, only the particles 
with polarization Jz = ±2 contribute to Xb2 decay into two transversely-polarized vector 
mesons. Similarly, one finds 



{0\xUb\Vb) 



N 

7r^.,(o,/i0. 

Ztc 



(44) 



The factorization of the matrix elements (l35l) implies that the decay rate also factorizes. 
For the decays of Xbo and Xb2 into two pseudoscalar mesons or two longitudinally-polarized 
vector mesons, we find 



r(x 



m 



bO 



AA) =-- 



D 



171: 



XbO 



'\l2 



mm 



XbO 



n ■ v' n ■ V 



2n 



|C(^)ni?;^„ (0,^)1 



dui du) ^ , 
— T 



n 2 



(45) 



UJ, fx, fx 



I. 3 



and 



r {xb2 AA) 



m 



D 



I — 4m r 



m 



Xb2 



DSN, 



/\|2 
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Xb2 



n ■ v' n ■ V 



2n 



duj duo 



n 2 



(46) 



— T [uj, UJ, fi, fi'; ^Pj) 0a(c^, f^')4>A{i^, /i') 



UJ UJ 



19 



where A = P,Vl- For the decay of Xb2 into two transversely-polarized vector mesons, one 
finds the decay rate by summing over the possible transverse polarizations: 



mi^^lnii, —Am? 



D\ "'x,2-'^"'D 3N, 



5 87rm^^2 27r 



Xb2 ' 

1 2 



F (/i)— — / —T[u;,uj,n,n; Pj) (pVj.{uj, h)(Pvt{uj, fJ^) 

Z Z I UJ UJ 



(47) 



In the case of r]}, decay into a pseudoscalar and a longitudinally-polarized vector meson, we 
find 



mi. ^ I ml. — 4m^ 



r ivb - PVl + c.c.) = ^ -^\C (/i) p \R,so, P ^ 



I n ■ V n ■ V 



— — T {uj,uj,n,n'; ^So) {(t)vi^{uj,^i')(j)p{uj,^i') - (j)v^{uj, ^i')(j)p{uj, ^i')) 

UJ UJ 



(48) 



Note that we are working in the limit mc oo, where the m£,*—m£) mass splitting vanishes. 

The factorized formulas Eqs. fl35l) and fH5l) - (jUj) are the main results of this paper. Each 
decay rate of (H5!) - (148|) depends on two calculable matching coefficients, C and T, and three 
non-perturbative, process-independent matrix elements, namely, two Z^-meson distribution 
amplitudes and the bottomonium wavefunction. In Sec. |V]we will provide a model-dependent 
estimate of the decay rates ( l45l) - ( 148|) and will discuss the phenomenological implications. 
We conclude this section by observing that all the non-perturbative matrix elements cancel 
out in the ratios r(xw PP)/^{Xb2 ^ PP) and T{xm ^ VLVL)/T{xb2 ^ V^Vl), since 
the spin symmetry of pNRQCD guarantees i?^^^^(0) = i?^^^(0), at leading order in EFTn. 
Neglecting the Xbo - Xb2 niass difference, we find, up to corrections of order w"^, 

T{xm ^ AA)/Tixb2 AA) = 1^ = 10 , (49) 

with A = P,Vl. 
B. Running 

The dependence of the matching coefficient T{uj, uj, /i, /x'; ^^'^^Lj) and of the operators in 
Eqs. ( |45|) - (j48l) on the scale /x' is driven by a RGE that can be obtained by renormalizing 
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the EFTn operators. The RGE for the EFTn operators, which also defines the anomalous 
dimension 7EFTn) is similar to Eq. ffTOl) . 

To calculate the anomalous dimension at one loop, we compute the divergent part of the 
diagrams in Figs. [H] and [71 As mentioned in Sec. [TTl the pNRQCD Lagrangian has the 
following structure, 

/ 73 rtusoft I 7- 

pNRQCD —lax -t-NRQCD + -^pot 5 

where the superscript usoft indicates that the gluons in the NRQCD Lagrangian are purely 
ultrasoft {■rrihw'^ , rrihw'^) , while Lpot contains four-fermions operators, which are non-local in 
space. 



At leading order in as{mbw) and r, is the Coulomb potential 



For the explicit form of higher-order potentials, see, for example, Refs. 12| 3l|]- Vertices 
from Lpot generate one-loop diagrams as the first diagram in Fig. [6l However, these diagrams 
do not give any contribution to the anomalous dimension at one loop. Indeed, the insertion 
of the Coulomb potential 1/r in Fig. [6] does not produce UV divergences. Insertions of the 
l/rrih potentials yield divergences but the coefficient of the l/rrif, potential is proportional to 
a'^{mhw), so it is not relevant if we are content with a NLL resummation. Insertions of 1/ml 
potentials give divergences proportional to subleading operators, which can be neglected. 
The second diagram in Fig. [6] yields a result completely analogous to the last term in Eq. 
f|T3|l . with the only difference of a color pre-factor, 

^-MpNRQCD = i^^^^N ■ (51) 

This divergence is completely cancelled by the 6-quark field renormalization constant Zb, 
and hence the pNRQCD diagrams in Fig. [6] do not contribute to the anomalous dimension 
at one loop. 
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FIG. 6: One-loop diagrams in pNRQCD. The first diagram contains insertions of quark-antiquark 
potentials. In the second diagram the gluon is ultrasoft. 





FIG. 7: One-loop diagrams in bHQET. There are three analogous diagrams for the other copy of 
bHQET. 

On the bHQET side, the third diagram in Fig. [7] is convergent, and hence it does not 
contribute to the anomalous dimension. The first two diagrams give 



iM 



bHQET,] 



(52) 



with 



a. 



A(u;, uj' , uj, uj') =—CfS{uj — uj') < 6 {uj — uj') 
2n 



2e^ 



+ 



1 



UJ — UJ' 



UJ 



d{uj- uj') +d{uj' -uj)e {uj) — 



(53) 



The diagrams for the bHQET copy in the n-direction give a result analogous to Eqs. (!52|) 
and fl53l) . with uj uj , uj' —>■ uj', and n ■ v' n ■ v. Extracting 7EFTn from the divergence is 
again standard, just as we did in the case of 7EFTr After adding to Eq. (153!) the bHQET 
field renormalization constants and for heavy and light quarks 

1 a. 



1 + -^c, 

e 271 



1 - -^C, 



we find 



7EFTii(t^,t^';t^,t^';/w') = 2'^p5 [uj - uj') 5{uj - uj') + 'yo{uJ,uj';uj,uj'; jj.') 



(54) 
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with 



7o(u;,u;';a;,a;';/i') 



a 



An 



- ^4CfS {uj - u') 

4:71 

- —4Cp6 (uj - u') 

4:71 



9{uj- iu') 

e{uj- Lj') 



-1 + ln 
1 

U — Lj' 
1 

UJ — Uj' 



UJ' 



In 



fi'n ■ V 



LU' 



UJ 



+ 9(uj-uj)9 (iu) — 



(55) 



UJ 



+ eiuj'-u)e{u)-i- 



1 



UJ \UJ' — UJ 



The term proportional to •yp in Eq- dSH) reproduces the running of POl) . 70 is 

responsible for the running of the D-meson distribution amplitudes and it agrees with the 
result found in Ref. [361]. Also, in Eq. (155|) the coefficient of In /i' is proportional to F 



cusp 



.a. 



Note that, since the bHQET Lagrangian is spin- independent, the anomalous dimension does 
not depend on the spin or on the polarization of the D meson in the final state, at leading 
order in the power counting. 

Using Eqs. fl50l) and fl54|) we find the following integro-differential RGE for the operator 
0{uj, LJ, /i'): 

d 



din /i 



-0{uJ,UJ,fi') 



duj' I cia;'7o(ijj, u;'; ijj, cj'; /i')C'(a;', tu', yu') , 



(56) 



where we have dropped both the subscripts A, B, and the superscript ^'^+^Lj, since 70 
does not depend on the quantum numbers of the initial or final state. Using the fact 
that the convolution of F'^(jj,')T{uj,uj, fx, fx'; ^'^+^Lj) and the operator {uj,uj, fx') is 

/i'- independent, we can write an equation for the coefficient. 



d 



din n' 



, UJ UJ 



duj' / duj' — —F (fi')T(uj',uj',fi,fi')'joi^\^]^\^]fJ'') 

' Uj' Uj' 



duj' / duj'F {fi')T{uj',uj',fi,fi')'yo{^,^']^,^']fJ'') 



(57) 



where the last line follows from the property of 70 at one loop. 



UJ UJ 



-— 7o(u;', uj; uj', u; /i') = 7^(0;, u'; uj, uj'; /i') , 
uj' uj' 

as can be explicitly verified from the expression in Eq. flS5]) . 



Eq. (1571) can be solved following the methods described in Ref. [36| . We discuss the details 
of the solution in App. [Al where we derive the analytic expressions for T{uj,u), fi, ^'; ^Pj) 
and T{uj,lD, fx, fx'; ^Sq), with the initial conditions at the scale fi'^ = rric expressed in Eq. (|33|) . 
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V. DECAY RATES AND PHENOMENOLOGY 



In Sec. IIV Al we gave the factorized expressions for the decay rates (H5l) - psj) : r(x5o, 2 
PP), T{xm,2 VlVl), T{xb2 ^t^t), and T{r], PVl + c.c.). In Sees. |IIlB]and[IVB]we 
exploited the RGEs ffT^ and (lFr|) to run the scales fi and /i', respectively, from the matching 
scales /i = 2mfe and fi' = rric to the natural scales that contribute to the matrix elements, 
/i = rric and yu' ~ 1 GeV, resumming in this way Sudakov logarithms of the ratios mc/2m;, 
and rric/ 1 GeV. 

We proceed now to estimate the decay rates fHSjl - fHHj) . In order to do so, we need to 
evaluate the following ingredients: the light-cone distribution amplitudes of the D meson and 
of the longitudinally- and transversely-polarized D* mesons, and the wavefunctions of the 
states rjf, and Xbj- In principle, these non-perturbative objects could be extracted from other 
rif,, Xb, and D-meson observables. In the case of the rjh, the value of the wavefunction at the 
origin can be obtained from a measurement of the inclusive hadronic width or of the decay 
rate for the electromagnetic process rji, —>■ 77, since they are both proportional to |i?,,j(0)p. 
Unfortunately, at the moment there are not sufficient data on rji, decays. Another way to 
proceed is to use the spin symmetry of the leading-order pNRQCD Hamiltonian, which 
implies Rti^{0) = Rj (0), and to extract the Upsilon wavefunction from r(T 



1.28 ±0.07 KeV 



371]. Using the leading-order expression for r(T e^e ) 



e^e ) = 
one finds 



|i?^(0)|2 = 6.92 ± 0.38 GeV^ where the error only includes the experimental uncertainty. 
The above value is in good agreement with the lattice evaluation by Bodwin, Sinclair, and 



Kim [39| and it falls within the range of values obtained with four different potential models, 
as listed in Ref. Q. 

1-^X60 2 ('-')l^ '"^^ obtained from the electromagnetic decay Xbo,2 — > 77. Unfortunately, 
such decay rates have not been measured yet. The values listed in Ref. 40] range from a 



minimum of \R' (0)]^ = 1.417 GeV^, obtained with the Buchmuller-Tye potential 4l|, to a 



XbJ ' 

maximum of ]i?^^^^(0)]^ = 2.067 GeV^, obtained with a Coulomb-plus-linear potential. The 
lattice value is roughly of the same size, \R'^^^{0)\'^ = 2.3 GeV^, with an uncertainty of about 
15% [39^]. We use this value in our estimate. 

For the pseudoscalar D-meson distribution amplitude we use two model functions widely 
adopted in the study of B physics. A first possible choice, suggested for example in Ref. 
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36l |. is a simple exponential decay: 

= IGeV) = ^^H^exp 



(58) 



Another form, suggested in Ref. 42|, is 



IGeV) 



UJ 



2{ao - 1) 



1 



CD2 



Incj 



7r^ 



(59) 



where u = u / jj! . The theta function in Eqs. (158|) and (159|) reflects the fact that the 
distribution amplitudes (j)A{^^, fJ''), with A = {P, Vl, Vt}, have support on u > [43]. 

The subscript indicates that these functional forms are valid in the D-meson rest frame, 
with a HQET velocity-label Vq = (1,0,0,0). With the definition we adopt in Eq. fl5Bl) . the 
distribution amplitude is not boost-invariant and in the bottomonium rest frame, in which 
the D meson has a velocity {n ■ v,n ■ v,0) ~ {mc/2mb, 2mb/mc, 0), it becomes 

, , 1 , 



■>p,o 



UJ 



(60) 



n ■ V \n ■ V 

as shown in App. [Bl Ad and ctd in Eqs. (158|) and (159|) are, respectively, the first inverse 
moment and the first logarithmic moment of the D-meson distribution amplitude in the 
D-meson rest frame, 



Ad^(/x') 



oo 



duj 



UJ 

^ duj_ 

UJ 



(t)pfl{uj,fl') 

In 



Furthermore we assume that the vector-meson distribution amplitudes (j)Y^{u) and (pvri^) 
have the same functional form as (f)p{uj), but with different parameters A^* , ct/j* and Xp)* , 



The D-meson distribution amplitude and its moments have not been intensively studied 
unlike, for example, the 5-meson distribution amplitude. Therefore, we invoke heavy-quark 
symmetry and use the moments of the i?-meson distribution amplitude in order to estimate 



the decay rate. However, the value o: 
sum rules, Braun et al. estimated 



" Ar is affected by a noticeable uncertainty. Using QCD 
m ^bLmI = IGeV) = 0.460 ± 0.110 GeV, where the 
uncertainty is about 25%. Other authors 4^ 45| 46|] give slightly different central values 
and comparable uncertainties, so that Xp falls in the range 0.350 GeV < Ab < 0.600 GeV. 
The first logarithmic moment ap is given in Ref. 42|, up = o'b{^J'' = 1 GeV) = 1.4 ±0.4. We 
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assume that the moments of the D*-meson distribution amphtudes fall in the same range 
as the moments of (j)p{uj). 

We evaluate numerically the convolution integrals in Eqs. fl45p - (HSj) . We choose the 
matching scales fib and fi'^ to be 2mh and rric respectively. Using the RGEs we run the 
matching coefficients down to the scales u = and /i' = 1 GeV. For the b and c quark 



masses we adopt the IS mass definition [47 1 



rubilS) = — = 4730.15 ± 0.13 MeV , 

m (61) 
nicilS) = — ^ = 1548.46 ± 0.01 MeV . 



The values of as at the relevant scales are [371] tts(2mb) = 0.178 ± 0.005, as{mc) = 0.340 ± 



0.020, and as{l GeV) ~ 0.5. With these choices, the value of g in Eq. (]A5|) is g{mc, 1 GeV) = 
-0.12 ±0.02. 

The decay rates T{xbj ^^) with A = {P, Vl, Vt}, (H5|) - (H7|) . depend on the masses of 



the XbJ and of the D mesons, whose most recent values are reported in Ref. [37(]. Since the 
effects due to the mass splitting of the Xbj and D multiplets are subleading in the EFT power 
counting, we use in the evaluation the average mass of the XbJ multiplet and the average 
mass of D and D* mesons: m^^j = 9898.87 ± 0.28 ± 0.31 MeV and = 1973.27 ± 0.18 
MeV. Therefore, the velocity of the D mesons in Xbj decay isn-v = n-v' = m^^j/mij = 5.02, 
with negligible error. The decay rate T{r]h PVl + c.c.) fHHj) depends on the mass of the 
r/b, which has been recently measured: m^^^ = 9388.9^2 3 ± 2.7 MeV 2^. The velocity of the 
D meson in the rji, decay is n ■ v = n ■ v' = mrj^/mD = 4.76, again with negligible error. 

The decay rate T{xbo PP) (l45l) . obtained with 0^^p and (j)^"^^^^ separately, is shown 
in Fig. [HI In order to see the impact of resumming Sudakov logarithms, we show for both 
distribution amplitudes the results with (z) the LL and NLL resummations and [ii) without 
any resummation at all. In the plots, we call the resummed results NLL-resummed, indi- 
cating that Sudakov logarithms are resummed up to NLL. For both distribution amplitudes 
the resummation does have a relevant effect on the decay rate. In the case of the 
resummation decreases the decay rate by a factor of 2 — 1.5 as Ad goes from the lowest to 
the highest value under consideration. In the case of 0^raun decay rate decreases too, for 
example, by a factor 1.5 when ar, = 1.4. In Fig. [9] we compare the decay rates obtained with 
the two distribution amplitudes. Over the range of A^i we are considering the two decay 
rates are in rough agreement with each other. 
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FIG. 8: T{xbo ~^ PP) as a function of \d, calculated with the distribution amplitudes (/j^^p (^e/i) 
and (^Bi'S'U'i {right). The dash dotted and solid lines denote the NLL-resummed decay rate. For 
comparison, the decay rate without resummation is also shown, denoted by dash double-dotted 
{left) and dashed {right) lines. For i^^"^^""^ we vary the parameter ao from ao = 1 (lower curve) 
to gd = 1-4 (middle curve) to (T£, = 1.8 (upper curve). 




FIG. 9: r(xf,o ~^ PP) a-s a function X^. The dash dotted line denotes the decay rate calculated 
with (f^'^^ , while the three solid lines with i;^^'''^"'^. For ^Braun vary the value of the parameter 
ud from = 1 (lower curve) to ao = 1-4 (middle curve) to ud = 1-8 (upper curve). 

Figs. [H]-[9]also describe the relation between the decay rate T{xbo — > VlVl) and Ad*. 
According to Eqs. fj46|) and P7|) . the processes Xb2 PP, Xb2 VlVl, and Xb2 VtVt 
show an analogous dependence on the first inverse moments of the light-cone distribution 
amplitudes, and they differ from Figs. [8] -[9] by constant pre-factors. Therefore, we do not 
show explicitly their plots. 

Qualitatively, Figs. [S] - El show a dramatic dependence of the decay rate on the inverse 
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moment A^. Using Eqs. (H5|) . (160|1 and (lA16p . one can show that when 0^'''^"° is used, the 
decay rate is proportional to A^'^, while it scales as X^~'^^ when we adopt (/)^^p, with g defined 
in Eq. (lASp . As a consequence, the decay rate drops by an order of magnitude when A^ goes 
from 0.350 GeV to 0.600 GeV. The particular sensitivity of exclusive bottomonium decays 
into two charmed mesons to the light-cone structure of the D meson — much stronger than 
usually observed in D- and S-decay observables — is due to the dependence of the amplitude 
on the product of two distributions (one for each meson) and to the non-trivial dependence 
of the matching coefficient T on the light-quark momentum labels uo and uo at tree level. On 
one hand, the strong dependence on a relatively poorly known quantity prevents us from 
predicting the decay rate T{xm DD). On the other hand, however, it suggests that, 
if the decay rate is measured, this channel could be used to better determine interesting 
properties of the D-meson distribution amplitude, such as Ad and (Jd- The viability of this 
suggestion relies on the control over the theoretical error attached to the curves in Fig. [8] 
and on the actual chances to observe the process Xb ~^ DD at current experiments. 

The uncertainty of the decay rate stems mainly from three sources. First, there are 
corrections coming from subleading EFT operators. In matching NRQCD -|- SCET onto 
pNRQCD + bHQET (Sec. IIV Al) . we neglected the subleading EFTn operators that are 
suppressed by powers of Aqcd/'"^c and w^, relative to the leading EFTn operators in Eqs. 
(IM]) and ([32]). In matching QCD onto NRQCD + SCET (Sec. UTTAl) . we kept only Jefti 
([6]) and neglected subleading EFTj operators, suppressed by powers of A and w"^. These 
subleading EFTi operators would match onto subleading EFTn operators, suppressed by 
powers of Aqcd/^^c and w"^. Using w"^ ~ 0.1 and Aqcd/"^c ~ 0.3, we find a conservative 
estimate for the non-perturbative corrections to be about 30%. 

Second, there are perturbative corrections to the matching coefficients C and T. Since 
as(2mf,) = 0.178, we expect a 20% correction from the one-loop contributions in matching 
QCD onto NRQCD -|- SCET. In the second matching step, similarly, the one-loop corrections 
to T{u,LJ, fi, fi'; "^^^^Lj) would be proportional to asirric) ~ 30%. We can get an idea of 
their relevance by estimating the dependence of the decay rate ( H5i) on the matching scales 
/ife and yU^. If the matching coefficients C and T and the anomalous dimensions 7efti and 
70(0;, cij', cj, cj'; /i') were known at all orders, the decay rate would be independent of the 
matching scales and /x^. However, since we only know the first terms in the perturbative 
expansions, the decay rate bears a residual renormalization-scale dependence, whose size is 
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FIG. 10: Left: Scale dependence of T{xbo PP) on the matching scale We vary from a 
central value = Irrij, (solid line) to a maximum of /x^, = 20 GeV (dashed line) and a minimum 
of /ifo = 5 GeV (dotted line). The dashed and dotted lines overlap almost perfectly. Right: Scale 
dependence of T(xbo PP) on the matching scale fj,'^. We varied /x^ from a central value of 
fj,'^ = rric (solid line) to a maximum of /x^ = 2.5 GeV (dashed line) and a minimum of /i^ = 1.2 GeV 
(dotted line). 

determined by the first neglected terms. 

In Fig. [To] we show the effect of varying fit between Arrih ~ 20 GeV and rrib ~ 5 GeV on 
the decay rate, using (j)^'^^^'^. The solid line represents the choice = '^f^b-, while the dashed 
and dotted lines, which overlap almost perfectly, correspond respectively to yU;, = 20 GeV 
and /ife = 5 GeV. The dependence on is mild, its effect being a variation of about 5%. We 
obtain analogous results for the decay rate computed with 0^^p, which are not shown here 
in order to avoid redundancy. 

On the other hand, even after the resummation, the decay rate strongly depends on /x^. 
We vary this scale between 1.2 GeV and 2.5 GeV and we observe an overall variation of 
about 50%. We expect the scale dependence to be compensated by the one-loop corrections 
to the matching coefficient T{u!,u, ^, fi'; ^Pj). This observation is reinforced by the fact 
that the numerical values of the running factors U{^b,f^) and V{^'^,^') (defined respectively 
in Eqs. (126|) and (]A6P ) at NLL accuracy are smaller than expected on the basis of naive 
counting of the logarithms. As a consequence, the next-to-leading-order corrections to the 
matching coefficient could be as large as the effect of the NLL resummation. In the light of 
Fig. dni the one-loop correction to T(u;, u;, /x, /x'; ^Pj) seems to be an important ingredient 
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for a reliable estimate of the decay rate. 

A third source of error comes from the unknown functional form of the D-meson distri- 
bution amplitude. For the study of the 5-meson shape function, an expansion in a complete 
set of orthonormal functions has recently been proposed and it has provided a systematic 
procedure to control the uncertainties due to the unknown functional form The same 
method should be generalized to the B- and D-meson distribution amplitudes, in order to 
reduce the model dependence of the decay rate. We leave such an analysis to future work. 

To summarize, the calculation of the one-loop matching coefficients and the inclusion of 
power corrections of order Aqcd/'^c appear to be necessary to provide a decay rate with 
an accuracy of 10%, that would make the decays Xbj D^D^, Xbj D^D^ competitive 
processes to improve the determination of and ctd, if the experimental decay rate is 
observed with comparable accuracy. 

We estimate the decay rate Tijiiy PVl + c.c.) (HHj) using 0^^p and 0^'''^™ for both 0p 
and 0Vl. In the limit nic — * oo, spin symmetry of the bHQET Lagrangian would imply the 
equality of the pseudoscalar and vector distribution amplitudes, (pp = and hence the 
vanishing of the decay rate T{rih PVl + c.c). Assuming spin-symmetry violations, the 
decay rate depends on (z) the two parameters Xd = {\d + Xdi) /'^ and 5 = {Xdi — ^d) / {^0 + 
Xdi), if (p^^^^ is used, and on (ii) three parameters \d, and \o-di — o"_d|, if cfp^^^^ is used. 

The two plots in the left column of Fig. [TTlshow the decay rate, computed with 0^^p, as 
a function of Xd with 6 adopting various values, and as a function of 6 with now being 
the parameter. In the right column, the decay rate computed with is shown. Since in 

this case the decay rate does not strongly depend on 6, we fix it at 5 = and we show the 
dependence of the decay rate on and |(Td£ —cr£)\. We "normalize" the difference between 
the first logarithmic moments by dividing them by cr = 2a d- 

The most striking feature of Fig. [11] is the huge sensitivity to the chosen functional form. 
Though a precise comparison is difficult, due to the dependence on different parameters, 
the decay rate increases by two orders of magnitude when we switch from 0^^p to ^B^'a-^'^. 
Once again, this effect hinders our ability to predict T{rib — > PVl + c.c.) but it opens up the 
interesting possibility to discriminate between different model distribution amplitudes. 

Using Eqs. dH]) and flA17D . we know that T{rib PVl + c.c.) goes like A^^"^^ when (p^'^P 
used or A^'^ when (f)^"^^^"^ used. Fig. [11] appears to confirm this strong dependence on Ad- 
The plots in the lower half of Fig. [TT] refiect the fact that the decay rate vanishes if one 
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FIG. 11: Left: r{r]b — > PVl + c.c.) as a function of Xd and 5, computed using exponential 
distribution amplitudes 4>^^ and 4>y^^- Right: r(r/b — > PV^: + c.c.) as a function of Xd and 
|(T£)£ — (T£)|/cj, computed with the Braun distribution amphtudes (p^^^^'^ and <py^^^^'^. 

assumes 4>p{uj) = 

We conclude this section with the determination of the branching ratios B{xbo PP) = 
T{xm PP)/T{xbo ^ hght hadrons) and ^ PVl + c.c.) = r(r/b ^ PVl + c.c.) /T{7]t ^ 
light hadrons). At leading order in pNRQCD, the only non-perturbative parameter involved 
in the inclusive decay width of the rji, is |-Rr,j^(0)p j4|, 

2Im/i(iSo)iV, 



T{r]i, light hadrons) 



\RvM\- 



(62) 



Therefore, B{r]b PVl + c.c.) does not depend on the quarkonium wavefunction and the 
only non-perturbative parameters in B{rih PVl + c.c.) are those describing the D-meson 
distribution amplitudes. 



For P-wave states, the inclusive decay rate was obtained in Refs. ^ 49|], where the 
contributions of the configurations in which the quark-antiquark pair is in a color-octet .S*- 
wave state were first recognized. In pNRQCD the inclusive decay rate is written as 
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FIG. 12: Branching ratios B{xbo ^ PP) iWt) and S(r/b 
computed using the distribution amphtude (^^''^u'l. 



PVl + c.c.) (right). The latter is 



5l| 



T{xbo ~^ light hadrons) 



1 3N, 



Xb 



lmf,CPo) + ^lmfsCS,)£ 



(63) 



where the color-octet matrix element has been expressed in terms of the heavy quarkonium 



wavefunction and of the gluonic correlator S, whose precise definition is given in Ref. [50|]. S 
is a universal parameter and is completely independent of any particular heavy quarkonium 
state under consideration. Its value has been obtained by fitting to existing charmonium 
data and, thanks to the universality, the same value can be used to predict properties of 
bottomonium decays. It is found in Ref. [sO^ S = 5.3 1^2 2- The matching coefficients in 



Eqs. ( l62l) and ( l63l) are known to one loop. For the updated value we refer to Ref. [52|] and 
references therein. For reference, the tree-level values of the coefficients are as follows ^: 



Im/iOS'o) = «^(2m,)7r 



2Nr 



Im/i(3Po) = 3a^(2^^)^ 



2Nr 



o 



(64) 

With the above parameters, we plot B{xbo PP) and B{rib PVl+c.c.) as a function of 
Xd and Xd, respectively, in Fig. [121 Over the range of A/) we are considering, B{xbo PP) 
varies between 4 ■ 10~^ and 4 ■ 10~^; it is approximately one or two orders of magnitude 
smaller than the branching ratios observed in Ref. 25(] for Xbj decays into light hadrons. 
B{rih — * PVl + c.c.) depends on the choice of the distribution amphtude. Choosing the 



parameterization 



tBraun 



it appears that, despite the suppression at \aD* — a^l = 0, 



B{rjb PVl + c.c.) assumes values comparable to B{xbo PP) even for a small deviation 
from the spin-symmetry limit. If is chosen, the branching ratio is suppressed over a 



32 



wide range of \aD' — ctdI- The branching ratio Bijii, — > PVl + c.c.) was first estimated in 



531]. The authors of [53| assumed that the exclusive decays into DD* dominate the inclusive 



decay into charm, r(?7;, — > PVl + c.c.) ~ ^ijlh cc + X). With this assumption, they 
estimated the branching ratio to be in the range 10^'^ < -^(r^f, PVl + c.c.) < 10^^. Our 
analysis shows that such an assumption does not appear to be justified in the range of A^j 
considered in Fig. [121 while it would be appropriate for smaller values of the first inverse 
moments, for example for Ad ~ 0.200 GeV if the distribution amplitudes are described by 

0Braun 

Our estimates indicate that observing the exclusive processes r^f, — DD* + c.c. and 
Xb DD would be extremely challenging. A preliminary analysis for Xb D^D^ 5^ 
suggests that the number of T(25') produced at BABAR allows for the measurement of a 
branching ratio B{xbo ~^ D^D^) ~ 10^^, which is two or three orders of magnitude bigger 
than the values in Fig. [121 An even bigger branching ratio would be required for the smaller 
T(2S') sample of CLEO. However, we stress once again the strong dependence of the decay 
rates on the values of the first inverse moments. In particular, our estimates rely on the 
relation Ad = Ab, which is valid in the limit of mb,mc oo; even small corrections to the 
heavy flavor symmetry, if they had the effect of shifting the value of Ad towards the range 
0.250 — 0.350 GeV, could considerably increase the branching ratios. 



VI. CONCLUSIONS 



In this paper we have analyzed the exclusive decays of the C-even bottomonia into a 
pair of charmed mesons. We approached the problem using a series of EFTs that lead to 
the factorization formulas for the decay rates (Eqs. fHSl) - fj48l) ). valid at leading order in 
the EFT power counting and at all orders in a^. We improved the perturbative results by 
resumming Sudakov logarithms of the ratios of the characteristic scales that are germane to 
the dynamics of the processes. 

The decay rates (l45l) - (148|) receive both perturbative and non-perturbative corrections. 
Perturbative corrections come from loop corrections to the matching coefficients C and 
T, which are respectively of order as{2'mb) ~ 0.2 and as{mc) ~ 0.3. The largest non- 
perturbative contribution could be as big as Aqcd/'^^c? which would amount approximately 
to a 30% correction. Therefore, corrections to the leading-order decay rates could be notice- 
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able, as the strong dependence of the decay rates on the renormahzation-scale /i^ suggests. 
However, the EFT approach shown in this paper allows for a systematic treatment of both 
perturbative corrections and power-suppressed operators, so that, if the experimental data 
require, it is possible to extend the present analysis beyond the leading order. 

For simplicity, we have focused in this paper on the decays of C-even bottomonia, in which 
cases the decays proceed via two intermediate gluons and both the matching coefficients 
C and T are non-trivial at tree level. The same EFT approach can be applied to the 
decays of C-odd states, in particular, to the decays T DD and T D*D*, with the 
complication that the matching coefficient T arises only at one-loop level. Moreover, the 
same EFT formalism developed in this paper can be applied to the study of the channels 
that have vanishing decay rates at leading order in the power counting, such as rjh —>■ D*D*, 
T DD* + c.c, and Xb2 DD* + c.c. Experimental data for the charmonium system 
show that, for the decays of charmonium into light hadrons, the expected suppression of the 
subleading twist processes is not seen. It is interesting to see whether such an effect appears 
in bottomonium decays into two charmed mesons, using the EFT approach of this paper to 
evaluate the power-suppressed decay rates. 

Finally, in Sec. |V] we used model distribution amplitudes to estimate the decay rates. 
The most evident, qualitative feature of the decay rates is the strong dependence on the 
parameters of the D-meson distribution amplitude. Even though this feature may prevent 
us from giving reliable estimates of the decay rates or of the branching ratios, it makes the 
channels analyzed here ideal candidates for the extraction of important D-meson parameters, 
when the branching ratios can be observed with sufficient accuracy. 
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APPENDIX A: SOLUTION OF THE RUNNING EQUATION IN PNRQCD + 
BHQET 



The RGE in Eq. can be solved by applying the methods discussed in Ref. [36| to 
find the evolution of the i?-meson distribution amplitude. We generalize this approach to 
the specific case discussed here, where two distribution amplitudes are present. Following 
Ref. [361], we define 

1 



ujT{uj, uj', as) 



Lange and Neubert 



with 



TT 



9{u- uj') 



UJ — UJ 



UJ 



+ e{u' -uj)9 (u) — 



UJ' \UJ' — UJ 



36| prove that 

duj'ujT{uj, uj', as)iuj'y = uj''J^{a, a,) 



(Al) 



s^F 
TT 



'l + a)+ij{l-a)+2jE] . 



ip is the digamma function and 'Je the Euler constant. Eq. (1A1I) is valid if —1 < Re a < 
1. Exploiting (lAip . a solution of the running equation Eq. (1571) with initial condition 
T{uj,uj,fiQ) = {uj/ fi'^y {uj/ft'o)^ at a certain scale /Xq is 

F\f^')T{uj,cu,ix')=F\li',)fiuj,ix',f^'„rj)f{cu,ix',ft'„0 , 



(A2) 



with 



UJ 



v-g 



f{uj, fi', /i'o, V) = (^] {n-vY exp U (/Iq, /x', t]) 
g = ^(^0,^') = / -r^rcusp(a) , 



7i(«,) = -2 



da 



Air 



cusp 



- + 7i(a)+J^(r/-^,a) 



(A3) 



The function f{uj, /i', /ig, has the same form as f{uj, /i', /Iq, r^) and is obtained by replacing 
UJ ^ UJ, 1] ^ ^, and n ■ v ^ n ■ v' in Eq. flA3p . The integrals over a can be performed 
explicitly using the beta function in Eq. ( 1281) . The result is 

" -9 



/(u;,/i',/io,r/)/(^,/i',/io,0 = ( — 



UJ \ f UJ 



^0 



[n ■ vn ■ V 



'\9 



exp [^(/^[i,^')] 



r(i - + ^)r(i + r^) r(i - ^ + ^)r(i + p 
r(i + - ^)r(i _ r(i + e - ^7)r(i - o 



(A4) 
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Where, at NLL, 



and 



T./ , M -rco) 27r Jr-l-rlnr / dulp /3i\ 1-r + lnr A s 



cusp 



+ ^(2-87i^)lnr 
Po 



(A6) 



with r = a;s(/i')/as(/iQ). Notice that in the running from /ig = to /i' = 1 GeV only three 
flavors are active, so in the expressions for /?o, A, and riusp we use n/ = 3. 

Eq. (lA4p is the solution for the initial condition T(c<j, cj, /Xq) = [u / ^'^f {uj/ fi'^)^. To solve 
the RGE for a generic initial condition, we express T as the Fourier transform with respect 

to InCj/yUg, 

T{uj, /io) = — ^ f drds exp T-ir In ^ exp (^-is In -^^ F[T] (r, s, /ig) 

where F\T] denotes the Fourier transform of T. From the solution (1A2I) - (1A4I) it follows that 

F^(/i')T(.;,c.,^')=^^ J ^ d^ds[yj i^^J (n ■ v n ■ v'Y F[T]{r, s, f,',) 

, r(l + ir + g)V{l - tr) T{1 + ts + ^)r(l - ts) 
exp [V{fio,fi)] 



r(l -ir- g)T{l + ir) T{1 - is- g)T{l + is) 

(AT) 

The Fourier transform of the matching coefficient in Eq. (\A7\i has to be understood in 

55l | . That is, we define the Fourier transform of T as the function 



the sense of distributions 
of r and s that satisfies 



(27r 



^ / drds F[T]{r,s, fx')LpAir, n')LpBis, n') = [ ——T{uj,uj,ij')(f)Aiuj,n')(f)Biuj,iJ,'), 
T^Y J Jo ^ ^ 

(A8) 

or, more precisely, F[T]{r, s, fi') is the linear functional that acts on the test functions ifA{f) 
and fB{s) according to 
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The function yj^ is the Fourier transform of the D-meson distribution amphtude, 



duj f uj 



(AlO) 



where the integral on the r.h.s. should converge in the ordinary sense because of the reg- 
ularity properties of the D-meson distribution amphtude. As in Sec. HVt the subscript A 
denotes the spin and polarization of the D meson. 

In the distribution sense, the Fourier transform of the coefficient 1/{uj + ui) is 



LU + LU 



(r, s, /i'o) = {27if^S{r + s + i) sech 



— r 
2^ 



S 



(All) 



where R = r + s, S = r — s, and the factor ^ comes from the Jacobian of the change of 
variables. The hyperbolic secant is defined as sech = 1/ cosh. Similarly, we find 



LU — U 
UJ + UJ 



(i?,5,/i'o) = ^(27r)25(i?) (^cosech 



TT 



-S + is 



+ cosech 



TT 



-5* — ie 



(A12) 



The S function in Eq. (lAlip has complex argument. The definition is analogous to the 



one in real space 
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{S{R + i),^{R)) = ^{-i) . (A13) 

Using Eqs. flAlip and flA12l) . we can perform the integral in Eq. flA7p . obtaining respectively 
T{u,uj, fi, fi'; ^Pj) and T{uj,lj, ^, ft'; ^Sq). In order to give an explicit example, we proceed 
using Eq. flAlip . Integrating the 6 function we are left with 



F\fx')T{uj, UJ- /i') = exp [V{^^',, /i')] - 



,,/2\ 1/2+S 
PO 



n ■ vn ■ V 



l\9 



dS exp 



S UJ 
-i—m — 

2 UJ 



sech 



jjQ \OJUJ ^ 

1 r(| + ^ + i5)r(| + ^-|5) 



(AM) 



i + ^2r(i-^-|5)r(i-^ + |5) • 

The integral flA14p can be done by contour. The integrand has poles along the imaginary 
axis. In S = ±i there is a double pole, coming from the coincidence of one pole of the 
hyperbolic secant and the singularities in 1/(1 + 5*^). The F functions in the numerator have 
poles respectively in 5* = ±z {2n + 3 + 2g) with n > 0, while the other poles of sech are in 
S = ±i{2n + 1), with n > 1. We close the contour in the upper half plane for a; > a; and in 
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the lower half plan for uo > uj, obtaining 



F\^^')T{u, /i') = F\^^',) exp [r(/x[„ ^')] 0{u - u)- 



1 / ^'qU ■ vn ■ v' 



l\ 3 



LOLO 



T{l + g)T{2 + g) 
T{l-g)T{-g) 



1 _ In - + - (7) - i){~g) + ^g)- ^(2 + g) 



1 r(l -n + ^)r(2 + n + ^) 



n=l 
oo 



cu/ ?2(n + 1) r(— 72 — 5f)r(l — (yf + n) 



n=l ^ ' 



1 



r(2 + n + 2^) 



[n + + n + ^) r(l + n)r(-n - 2^) 



+ (cj ^ cj) 



(A15) 



with csc(5f7r) = 1/ sin(5f7r) and i\} is the digamma function. More compactly, we can express 
Eq. (lAlSP using the hypergeometric functions 4F3 and 3F2, 



exp ^{[i^, /i')] 



yu/n ■ vn ■ V ^ 

uouo 



d{uo-uj) \ T{l + g)T{2 + g) 
u \ T{l-g)T{-g) 

lu r(^ + 2)r(^ + 3) 

2a; T{l-g)T{2-g) 

. . .r(2 + 2^)2 



1 _ In - + - (7) - ^(-^7) + ^(1 + - ^(2 + (?) 



+ 



4F3 (1, 1, ^7 + 2, (7 + 3; 3, 1 - 2 - g- 



4 cos(^7r) ^^" ' 3i^2 (g + l,2g + 2, 2^7 + 3; 2, (7 + 3; --) \ + {u ^ uj) 
g + 2 V cu/ J 



(A16) 

where we have introduced the constants that appear in the initial condition in Eq. (1331) . In 
the same way, we obtain 



T?2/ l\rp/ - I In \ Tp2/ f\^F 47rQ;s(/i') 



2m TTlb 



exp [V{fi'^,fi')] 9{uj-uj) 



lj,\rn ■ vn ■ V ^ 



LOLO 



T{l + g)T{2 + g)uj 



- (^) 4 cos(^7r)r(l + 2^)r(2^ + 2) ^Fi {2g + 2,2g + 1; 2; | 



[LO LO) 



(A17) 



In Eqs. (IA16P and (lAlTp we renamed the initial scale /Iq = fi'^ to denote its connection to 
the scale m^. Setting fi' = fi'^ or, equivalently, g = 0, it can be explicitly verified that the 
solutions Eqs. flA16P and flA17l) satisfy the initial conditions Eq. fl33|) . 
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APPENDIX B: BOOST TRANSFORMATION OF THE L>-MESON DISTRIBU- 
TION AMPLITUDE 

We derive in this Appendix the relation between the distribution amphtudes in the D- 
meson and in the bottomonium rest frames, as given in Eq. fIBUl) . In the D-meson rest frame, 
characterized by the velocity label fo = (1,0,0,0), the local heavy-light matrix element is 
defined as 

mlM^l^KmD)., = -^F{^^')^ . (Bl) 

The matrix element of the heavy- and light-quark fields at a light-like separation Zq = 
n ■ defines the light-cone distribution (poi^ ■ zo,fi') in coordinate space: 

mlin . zo) ^7,n^M\DU = -^F{^^')^Mn ■ ^o,/^') • (B2) 

Eqs. (IBI]) and imply 0o(O,/i') = 1. In the definitions and the subscript 
is used to denote quantities in the D-meson rest frame. This convention is used in the rest 
of this Appendix. In the bottomonium rest frame, where the velocity label in light-cone 
coordinates is f = {n ■ v,n ■ v,0) and the light-like separation is 2;'^ = n ■ zn^/2, we define 

mlw^i^Kmo)^ = -tFif,')^ (B3) 

and 

miin ■ z) ^l,n^M\D). = -zF(/x')^0(n ■ z,^,') . (B4) 

Suppose that A is some standardized boost that takes the D meson from f , its velocity 
in the bottomonium rest frame, to rest. It is straightforward to find the relations between 
the D-meson momenta in the two frames: 

n ■ po = n ■ vn ■ p and n ■ po = n ■ vn ■ p . 
There is a similar relation for the light-cone coordinates, 

n ■ Zq = n ■ vn ■ z . 

With f/(A), the unitary operator that implements the boost A, one can write 

U{A)\D), = \D),, . 
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We choose A such that, for the Dirac fields, 

UiA)dix)U-\A)=A-l,^\Ax) and f/(A)/i^(a;) [/-^(A) = A^/^ /i^(Ax) 

where 

Ai/2 = cosh - + smh - , 

with a related to f by e" = n ■ f and e~" = n ■ v. 

Now we can write the matrix element in Eq. (]B3[) as 

= {0\U-\A) (u{A)am-'W) {U{Amo)U-\A)) U{A)\D% 
= m\0)Ay,KA-l,h%0)\DU = n-v {0\^\o)Kh^{0)\D)^^ 



(B5) 



-IP [jj, ) n ■ Vq = —tr [fi 



2 ^' ' 2 

where, in the last step, we have used n ■ vq = 1. Eq. fIBSp is thus in agreement with the 
definition in Eq. fIBSp . Applying the same reasoning to Eq. flB4p . one finds 



VO 



(B6) 

Comparing Eq. flB6p with flB4p . we see that (j){n ■ z, fi') = (^^{n ■ vn ■ z^ fi'). Note that in the 
bottomonium rest frame the normalization condition for the distribution amplitude is also 
0(0, /i') = 1. 

In the main text of this paper we have used the D-meson distribution amplitudes in 
momentum space, 

M^o, f^') = ^ jdn- zo e'"""-^Vo(n ■ zo, /i') , 
0(cu,/x') = ^ [ dn- ze^^"->(n ■ z,fi') . 



27r „ 

Using Eq. (IB6p . we can relate the two distributions: 

0(cu, IJ,') = — [ dn-z e*^"-"0(n ■ 2, /i') = — I dn ■ z e*'^"-"0o(^ -vn-z^fj.') 

271 I 271 



dn ■ z e^^-^^'^cpoin ■ z, /i') = - — 0o (- — ; /^'l 

n ■ V \n ■ V / 



271 n ■ V i n ■ V \n ■ V 
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as stated in Eq. (160|) . The D-meson light-cone distribution is normahzed to 1 in both 
frames, 

duJo(f)o{uJo, fi') = / du(f){u,^') = 1, 



as can be easily proved using 0o(O, /u') = 0(0, fi') = 1. 
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